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ABSTRACT 

We investigate two and three-dimensional shell-structured-inflatable froths, which 
can be constructed by a recursion procedure adding successive layers of cells around a 
germ cell. We prove that any froth can be reduced into a system of concentric shells. 
There is only a restricted set of local configurations for which the recursive inflation 
transformation is not applicable. These configurations are inclusions between succes- 
sive layers and can be treated as vertices and edges decorations of a shell-structure- 
inflatable skeleton. The recursion procedure is described by a logistic map, which 
provides a natural classification into Euclidean, hyperbolic and elliptic froths. Froths 
tiling manifolds with different curvature can be classified simply by distinguishing be- 
tween those with a bounded or unbounded number of elements per shell, without any 
a-priori knowledge on their curvature. A new result, associated with maximal ori- 
entational entropy, is obtained on topological properties of natural cellular systems. 
The topological characteristics of all experimentally known tetrahedrally close- packed 
structures are retrieved. 

I. INTRODUCTION 

A froth is a (topologically stable) division of space by cells, which are convex 
polytopes (polygons in 2D, polyhedra in 3D) of various shapes and sizes. These geo- 
metrical systems have attracted much attention in recent years, both theoretically and 
experimentally [1], [2]. The aim in this work is to study a specific class of froths, 
namely those which are reducible to a set of concentric shells. These particular froths 
are structured as if constructed in the following way. In a first stage, cells are added 
to a germ cell, forming around it a first layer whose external surface constitutes the 
second shell. In a second stage, cells are added to the first shell so as to form a second 
layer of cells encircling the first one, and so on. We emphasize that the words "germ" 
and "stage" are purely pictorial and do not imply any particular mode of growth since 
any cell of a generic shell-structured froth may play the role of its germ cell. Such a 
froth is called shell-structured-inflatable from now on. 
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A definition of a shell-structured-inflatable froth requires the notion of topological 

distance between cells. The topological distance t between two cells A and B is defined 
as the smallest number of edges crossed by a path connecting A and B. The germ cell 
is therefore at distance t — 0. A shell (t) is defined as the interface between two sets of 
cells distant by t and t+1 from the germ cell. A 2D froth is shell- structured-inflatable 
if it satisfies the following two conditions : 

1) For any set cells, equidistant to the germ cell, there exists a closed non self- inter- 
secting path which goes only through these cells and connects all of them. 

2) Any cell at distance t from the germ cell is the neighbour of at least one cell at 
distance t + 

Two consecutive shells (t) and of a shell-structured- infiatable froth are connected 

through a set of disjoint edges with one vertex on shell [t) and the other on shell (t + 1). 
These two shells are closed loops of edges delimiting the layer {t+ 1) of cells which are 
at distance t + 1 from the germ cell. Shell (t) divides the froth into an internal froth, 
constituted of cells at distances r < t, and an external froth, with cells at distances 
r > t. The extension to a 3D shell-structured-infiatable froths is straightforward and 
is given in Appendix B.2. 

In this paper we prove that the 2D and 3D shell-structured-inflatable froths are 
constructed according to a recursion procedure which is the logistic map [3], well- 
known in the theory of dynamical systems. The logistic map provides a natural clas- 
siflcation of these froths according to the behaviour of the number of edges per shell 
as the topological distance t increases. 

Any given froth is not necessarily shell-structured-inflatable. However, it has to be 
noted that a froth can always be decomposed into shells with respect to an arbitrarily 
chosen germ cell. In this decomposition, each cell of the layer (t) belongs to one of two 
categories. The cells of the flrst category, individually, have neighbours in both layers 
(t — 1) and (t + 1) and, collectively, are building up a complete ring around the chosen 
germ cell. The set of all these rings constitutes the "skeleton" of the shell-structure. 
The cells of the second category have neighbours in only one of the two layers (t — 1) 
or (t + 1). These cells can be considered as local topological defects included between 
the rings of the "skeleton" of the shell-structure. The "skeleton" is itself a space-fllling 
froth which is shell-structured-inflatable. The recursion procedure that we are studying 
applies to such a structure. 

The plan of this paper is the following. In Section 2, we derive the recursion pro- 
cedure associated with 2D shell-structured- inflatable froths and show that it can be 
written as the logistic map. The resulting classiflcation into Euclidean, hyperbolic and 
elliptic froths is discussed. In Section 3, it is shown that the recursion procedure in 3D 
is again described by a logistic map. The curvature of the embedding space is classi- 
fied as for the 2D froths. Section 4 gives examples of space-filling cellular structures 
which fit into the classification of 3D shell-structured-infiatable froths provided by the 
logistic map. In Section 5, a new bound on topological properties of natural cellular 
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structures is obtained. The topological properties of all experimentally known tetra- 
hcdrally close-packed (t.c.p.) structures are retrieved under the hypothesis of shell- 
reducibility. A conclusion emphasizes the main results of the paper. In Appendix A, 
the recursion procedure is generalized to 2D shell-structured- inflatable networks with 
vertex coordination larger than 3. Local topological defects in 2 and 3D shell- reducible 
but not inflatable froths are considered in Appendix B. Random 3D Euclidean froths 
are constructed from 2D random shell networks in Appendix C. 

II. RECURSION PROCEDURE FOR 2D FROTHS 

This section is concerned with two dimensional shell-structured-inflatable froths. 
The recursion procedure is derived here for froths and it is extended to networks with 
vertex coordination larger than 3 in Appendix A. Fig.(l) shows an example of a froth 
with the various shells indicated by bold lines and labelled by the index t (the shell 
t = corresponding to the boundary of the germ cell). Let K^('L) be the number 
of vertices going out from shell (t) towards shell {t + 1) (towards shell {t — 1)). Let 
F^^^ be the number of cells in the layer between shells (t) and {t + 1). If (n) is the 
average number of edges per cell in layer (t), the edges in this layer are accounted for, 
as follows 

(n)F(^) = yi*) + 2vi*^ + 21^*+') + Vi'+'^ . (2.1) 

In the right-hand-side of this equation, the quantity v!i^ + V^^ is the total number 

of vertices constituting shell (t), the quantity vl*^^"* -|- V^~^^^ is the total number of 

vertices constituting shell {t + 1) whereas the quantity v!l~^^^ + V^^ gives the number 
of vertices (counted twice) bounding the edges separating the cells comprised between 

shells (t) and {t + 1). Since = and F^*) = one has the recursion 
equation 

{n)vi*^ = 4F|*) + + Vi*^ . (2.2) 
The matrix form of this recursion equation is 

-.■)(?) ■ 

with recursion parameter s = (n) — 4. Eq.(2.3) generates recursively the whole froth 
from the germ cell. In general, the quantity (n) changes from one layer to the next. 
Hence the recursion parameter should depend on the distance t. However, the value 
of (n) associated to a layer of cells at a distance t from the germ cell must, as t — > oo, 
converge to the average value for any cell in the froth. Moreover, since the choice of 
the germ cell is completely arbitrary, the quantity (n) associated with layer (t) is an 
average. Consequently, the recursion parameter can be taken as an effective quantity 
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which is independent of t, and the quantity (n) is then the average number of edges per 

cell in the froth. The initial conditions in Eq.(2.3) are then v!^^ = and = (n). 

The recursion procedure described by Eq.(2.3) appears also in other instances, 
such as in the computation by decimation of the electronic energy spectrum in the 
ID tight- binding model [4]. In this case, the variables V^^^ are replaced by the 
components of the electronic wave-functions in the basis of the site states, and the 
recursion parameter s is the (dimensionless) energy of the electron. 

Eq.(2.3) gives an immediate link between the shell-structured-infiatable froths 
and the logistic map. Indeed, from the relations sV^^ = V^~^^^ + ^\ sV^~^^^ = 

s,vl'^ = vl'+'^ + vt'^ , (2.4) 

with 

, (2.5) 
and a similar relation for the V_'s. Iterating j times, one obtains 

s,vi'^ = vl'+''^ + vf^ , (2.6) 

with 

= s| - 2 , (2.7) 

and So — s. Eq.(2.7) is the trace map of the recursion matrix in Eq.(2.3). It is a 
logistic map [3], with two (unstable) fixed points s* = 2 and s* = —1. The logistic map 
decomposes the axis of values of the recursion parameter s into two different regions. 
Any point in the region |s| > 2 is sent towards infinity by the successive iterations 
of the logistic map. By contrast, if \s\ < 2, successive iterations of the logistic map 
remain all within this interval. The existence of these two intervals classifies all 2D 
shell-structured-inflatable froths. This classiflcation, corresponds to the curvature of 
the manifold which the froth tiles. The space is elliptic for \s\ < 2, hyperbolic for 
\s\ > 2 and Euclidean for the fixed point s = s* = 2. The map relate successive 
numbers V^*^ of vertices per shell. Iterations of Eq.(2.3) generate trajectories in the 

plane (t, V+) starting from the initial points V^T^'^ — = and = (n). 
When |s| < 2, the trajectories are given by the equation 

sin 

with cos((/?) = s/2. Eq.(2.8) shows that all trajectories are finite and end at the point 
= 0, with T = ^ — 1- Moreover, the values of V^^ are bounded by the quantity 

/sm{ip). These finite and bounded trajectories are describing the iterative tiling of 
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compact manifolds with positive curvature. Indeed, consider a froth tihng the surface 
of a sphere. Suppose that the north pole of the sphere is located in the germ cell; 
the successive shells are the parallels on the sphere. The number of vertices per shell 
increases between the north pole and the equator, then decreases from the equator 
to the south pole where the tiling ends. This is precisely the behaviour described by 
Eq.(2.8). The quantity T + 1 = ^ is the topological distance between both poles. Here 
are a few examples of regular froths with |s|<2. Tos = — 1 corresponds to a froth 
made with four triangles, i.e. the surface of a tetrahedron. The recursion parameter 
s = corresponds to a froth made with six squares, i.e. the surface of a cube. Finally, 
s = 1 is associated with a froth which is the surface of a dodecahedron. 
In the case \s\ > 2, the solution of Eq.(2.3) is 

yit)^yiO)^^M^it + l)) (2.9) 

sinh ip 

with cosh(</?) = s/2. Eq.(2.9) shows that, contrary to the previous case, the values of 
V^"^ increases exponentially with t. All trajectories are now infinite and unbounded 
in the plane (t,V+). They are therefore describing the iterative tiling of non-compact 
manifolds with negative curvature. 

At the fixed point s = s* = 2, Eq.(2.3) has the solution 

= (t + 1)K[°) , (2.10) 

The values of V^'^ have again no upper bound, but here they are increasing linearly 
with t as expected for the Euclidean plane by simple geometrical considerations. The 
fixed point s* = 2 describes shell-structured-inflatable froths covering the Euclidean 
plane with cells with 6 edges on average. An example of such froths is the hexagonal 
tiling. 

We have shown that the logistic map provides, in a natural way, the topological 
classification of tilings of manifolds without any a-priori knowledge of their Gaussian 
curvature. In 2D this classification by the logistic map is identical to that provided by 
the combination of the Gauss-Bonnet theorem [5] and Euler's equation: 

Kda=^{6-{n))F=^{2-s)F , (2.11) 

here, k is the Gaussian curvature and it is integrated over the whole manifold. F is the 
total number of cells in the manifold. The tiled manifold is hyperbolic, Euclidean or 
elliptic according when the integrated curvature is negative, zero or positive, i.e. when 
the recursion parameter s is larger, equal to or smaller than 2. However, the logistic 
map is also applicable in 3D where there is no Gauss-Bonnet theorem and the Euler 
equation is homogeneous [6]. 
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III. RECURSION PROCEDURE FOR 3D FROTHS 



This Section extends the analysis of the previous one to 3D sheU-structured- in- 
flatable froths. The froth has V vertices, E edges, F faces and C polyhedras. Every 
shell of the 3D froth is built up from two superposed different two- dimensional froths, 
and looks like a corrugated sphere. This is the same as in 2D, where a shell can be 
regarded as the superposition of two ID froths, one whose vertices are connected to 
the "incoming" edges from shell {t — 1) to shell {t), and the other, whose vertices are 
connected to the "outgoing" vertices pointing from shell (t) towards shell (t+l). Simi- 
larly, every spherical shell (t) of the 3D froth is built up of the superposition of two 2D 
froths, one whose edges are connected to the "incoming" faces of layer (t — 1), and the 
other whose edges are connected to the "outgoing" faces of layer (t). Let K^('L) and 

be the numbers of vertices and edges of shell (t), bounding the cells layer (t) 
between shells {t) and {t+ 1) (layer {t — 1) between shells (t) and {t — 1), respectively), 
which are making the "outgoing" ("incoming") froth. Let F^^^_^ be the number of 
faces of such froths. Both froths are characterized by the identities 



(Euler's formula) and 



V^t) =2 , (3-1) 

3V|*)_) = 2E^+l) (3.2) 

(since in both 2D froths, any vertex is connected by three edges and any edge is 
bounded by two vertices). 

One has the following relations between two successive shells 

yit+l) ^ y{t) 

Shell {t) is a spherical surface tiled by a network with F^p faces. One has 

Fii+'^ = - 2£;J) - Fii^ . (3.4) 

In this equation, (/) is the average number of faces per cell in the layer (t). 

Since the whole shell (t) is a polyhedron, both Euler's formula and the incidence 
relations are applicable between the number of edges the number of vertices V^*'' 
and the number of faces F^^ of the shell- polyhedron, namely, 

l4*^-i?^J)+FW=2 (3.5) 
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and 

(n)iv4') = 2E(J) . (3.6) 

Here, (n) n is the average number of edges per face of the sheU- polyhedron. Since it 
is an eUiptic tiling with vertex coordination > 3, then (n) n < 6. The shell-network is 
the superposition of two 2D froths, it has therefore 3-connected vertices corre- 
sponding to the "outgoing" ("incoming") froth, and also 4- connected vertices l^x at 
the intersections between edges of the two 2D froths. The three types of vertices are 
represented in Fig. (2). Fig. (3) shows the shell- network in the particular case of the 
"Kelvin froth" [7] , [8] , and indicates a 3-connected vertex and a 4-connected vertex. 

The total number of vertices V^*'' on shell (t) is the sum of all 3- and 4- connected 
vertices, i.e. 

y^*) = vl*^ + V^*^ + v!f^ . (3.7) 
The total number of edges E"^"* on shell (t) satisfies the equation 

2£;JJ^ = 3F|*^ + 3Vi*^ + 4F^*^ . (3.8) 
Using Equs.(3.5), (3.6) and (3.8), one obtains 

2 + Vl - (3y[*) + 3Fi*^ + 41/^*^) . (3.9) 

2 V {n)N ^ 



Combining Equs.(3.7) and (3.9), it is possible to express the variable V^^ in terms of 
the variables V^^ and V^^ alone as 

which, with Eq.(3.6) and (3.8), yields 

it) 8-(W*^ + y^*^) , , 

Putting Eq.(3.11) into Eq.(3.2), we obtain, with the help of Eqs.(3.1), (3.3) and (3.4), 
the following relation 

^^^'^ = ^(((/)-6)((n)Ar-4)-4)yf -F5 + 2(8 + (/)((n)A,-4)) . (3.12) 
Finally, by shifting the variables Vl^.(_) as 

8-((/)-6)((n)^-4) ' ^'-^'^ 



one obtains the recursion equation 



sVi'^ = vi'+'^ + vi'-'^ , (3.14) 

with the recursion parameter 

«=^(((/)-6)((n)^-4)-4) . (3.15) 
This recursion equation has the same matrix form as in the 2D case 




(3.16) 



As in 2D (/) and {n)N can be supposed to be independent of the distance t. The 
variation of the 3D recursion parameter s (the trace of the transfer matrix Eq.(3.16) ) 
is described by the logistic map (2.7), as in the 2D case. Consequently, the classification 
of 3D shell-structured-inflatable-froths is the same as in 2D. 

Elliptic shell-structurcd-inflatable froths are associated with |s| < 2. They are 
tiling iteratively 3D compact manifolds with positive curvature. Indeed, the corre- 
sponding solution of Eq.(3.16) is finite and bounded in the {t, V) plane 

=^sin(^t + S) + 2(^±i4^Mi^) , (3.17) 

with cos(<^) = s/2. The coefficients A and B can be deduced from the initial conditions 
= 2((/) - 2) and = 0. 

Hyperbolic shell-structured-inflatable froths are associated with |s| > 2. They 
are tiling iteratively 3D non-compact manifolds with negative curvature. Indeed, the 
corresponding solution of Eq.(3.16) is unbounded in the (t, F) plane 

= A sinh(^t + i?) + 2 ( ' + ^^^^"^^ - ) , (3.18) 

with cosh((^) = s/2. As previously, the coeflBcients A and B can be determined from 
the initial conditions. 

For s = — 2 the solution reads 

v!i^ = {-l)\At^B) + ^-±^^^MjiJl^ . (3.19) 
with A and B deducible from the initial conditions. 
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The solution of Eq.(3.16) associated to the fixed point s = s* = 2 is 

= (t + 1) (vl'^ +t(8+ {f){{n)N - 4)) j . (3.20) 

The quadratic dependence in t is the one expected from simple geometrical reasoning 
for a tiling of the 3D Euclidean space. 

As in 2D, the logistic map gives a natural description of tilings of three dimen- 
sional manifolds without the need of any a-priori information on their curvature. Con- 
sequently, the logistic map is able to characterize curved manifolds even when the 
Gauss-Bonnet formula is not applicable [9], [10]. The generation of tilings of curved 
manifold by the recursion procedure has therefore a wider applicability than the Gauss- 
Bonnet formula. 



IV. EXAMPLES OF 3D SHELL-STRUCTURED 
INFLATABLE FROTHS 



In order to illustrate the previous considerations, we give some known examples of 
3D froths and show that they fit our classification. All are monotiled (i.e. constituted 
of topologically identical cells), apart from the last example. 

The only regular elliptic froths in 3D are {3, 3, 3} (packing of tetrahedra), {4, 3, 3} 
(packing of cubes) and {5,3,3} (packing of dodecahedra) [11]. They correspond to 
s = — 1, s = — 2 and s = 1 respectively. Note that the case s = does not correspond to 
any regular froth. Indeed, the only solution s = of Eq.(3.15) with (/) and (n)jv < 6 
being both integers is (/) = 10, {n)^ = 5; which is not regular. 

Consider Eq.(3.15) in the Euclidean case (i.e. s — 2). This equation gives a 
relationship between the average number of neighbours per cell ((/)) in the 3D froth 
and the average number of edges per cell {{n)^) in the 2D spherical shell-network 

(/) = 6+^y^ • (4-1) 
(n)iv - 4 

This equation gives the condition for Euclidean space-filling by a shell-structurcd- 
inflatable froth. Note that, from Eq.(4.1), the minimal number of faces per cell of such 
a froth is 10, since (n) n < 6. It is known that the minimal number of neighbours 
per cell is 8 for an Euclidean froth. Thus an Euclidean froth with 8 < (n) tv < 10 
necessarily contains local topological defects of the kind shown in Fig.(ll). 

Recall that the shell-network is the superposition of two elliptic 2D froths, the "in- 
coming" and the "outgoing" froths. The pattern of edges constituting the shell-network 
sets the value of {n)N- Therefore, Eq.(4.1) allows us to construct systematically 3D 
Euclidean shell- structured-infiatable froths starting from the 2D shell-networks. 
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The simplest 2D froth is the hexagonal lattice. The examples displayed in Figs. (4), 
(5) and (6) illustrate the construction of ordered, monotiled 3D froths from a shell- 
network generated by different superpositions (cf. Figs. (4a), (5a) and (6a)) of two 
hexagonal lattices. Fig. (4b, c) shows two 3D unit cells constructed from the network 
(4a) (see also [12]). The cell (4b) is topologically equivalent to Kelvin's a- tetrakaidec- 
ahedron [7] [8] (it builds up the Kelvin froth shown in Fig. (3)), and the cell (4c), to its 
twisted variant [13]. Both structures have (/) = 14 and (n)jv = 5, Eq.(3.15) gives 
s = 2. They are indeed Euclidean space-fillers. 

Fig. (5a) shows part of a shell- network with 5-sided faces, generated by the super- 
position of two "squeezed" hexagonal lattices (see also [12]). Figs. (5b) and (5c) show 
two 3D unit cells constructed from the network (5a). These cells have again (/) = 14. 
The unit cell (5b) is topologically equivalent to the /3-tetrakaidecahedron (the Williams 
cell [14]). It has (n)Ar = 5, and is an Euclidean space-filler according to Eq.(3.15). 

The unit cell of Fig. (5c) is topologically equivalent to the 14-sided cell (the Gold- 
berg cell [15]) which occurs, among others, in clathrates [8], in t.c.p. structures [16] 
[17] and in the minimal froth of Weaire and Phelan [18]. Space can be filled layer by 
layers of Goldberg cells only. The layers (Fig. (5) ) are Euclidean and the network (5a) 
is the same as that of the Williams space-filler. However, successive layers are more 
and more distorted [19], as shown in Fig.(5d). This distortion, which stretches the 
network in one direction and compresses it in the other, strongly suggests that we are 
filling hyperbolic 3D space with a stack of Euclidean layers. It is possible to prove this 
contention by filling space shell by shell instead of layer by layer. When doing so, one 
finds that most of the shell-network is composed of pentagons (12 out of 14 in each 3D 
cell), but a finite density of hexagons (2 out of 14 in each 3D cell) is needed in order 
to close a shell. Thus, (n)7v > 5 which, according to Eq.(3.15), implies s > 2. Hence, 
the 3D manifold tiled by Goldberg cells is hyperbolic. 

With another intersection of the two "squeezed" hexagonal lattices, one generates 
the shell-network shown in Fig. (6a). The corresponding 3D unit cell (6b), has (/) = 16 
(8 quadrilaterals, 6 hexagons and 2 octagons) and (n)iv = 4.8. As far as we know, this 
unit cell is a new monotile Euclidean space-filler. 

Fig. (7) shows an example of an Euclidean shell-structured-inflatablc froth made 
of two different cells [8]. The shell-network (7a) has also two different tiles. The 
associated 3D unit cell (7b) has (/) = 12. 

One can show in general that any Euclidean shell-structured-inflatable froth made 
with topologically identical cells can be constructed from a shell-network generated by 
superposition of two hexagonal lattices. 

The construction of 3D disordered froths from 2D disordered shell-networks is 
discussed in Appendix C. 

Although a construction of 3D froths layer by layer has been given in [12], it must 
be emphasized that our approach, combining spherical shells with the logistic map, 
is more general and provides a unifying way to deal with 3D space-filling structures, 
whether regular or not, whatever the curvature of the manifold which they are tiling. 
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V. BOUNDS ON TOPOLOGICAL PROPERTIES OF NATURAL 
CELLULAR SYSTEMS AND T.C.P. STRUCTURES 

The average number ((n)) of edges per face of a 3D froth is in general different 
from the average number of edges per face in the shell- network ((n)Ar). For example, 
the froths in Figs. (4) and (5b) have {n)N = 5 and (n) = 5.14, the froth in Fig. (6) has 
(n)jv = 4.8 and (n) = 5.25 and the froth in Fig. (7) has {n)]\f = 5.33 and (n) = 5. 

The value of (n) is related to the average number of faces per 3D cell by 



It is interesting to study the competition between Eq.(5.1) and the Euclidean space- 
filling condition given by Eq.(4.1). These two relations (/)((?^)7v) (labelled "space- 
fiUing") and (/)((n)) (labelled "3D cell") are plotted in Fig.(8). They meet at the 
point ((n)*, (/)*) given by 



It is only when the equality (n) = {n)N = {n)* (which corresponds to (/)* = 13.29...) 
is satisfied that an arbitrary cell has the freedom to adhere to a preexisting shell by any 
subset of it faces, without adjustment. This freedom grants therefore a larger number of 
possibilities for building up a froth and it maximizes the orientational entropy per cell. 
Indeed, Eq.(5.1) is a constraint on any single 3D cell, whereas Eq.(4.1) is a constraint 
on the set of 3D cells in a layer. When (n) = {n)N = {n)*, one of the two constraints 
is automatically satisfied by the other and the orientational entropy is increased [20]. 

Note that the value (/)* = 13.29... falls within the range of several already known 
bounds. It is consistent with the values 13.2 and 13.33... resulting from the decurving 
of the dodecahedral packing with 14- and 18-sided cells or 14^ and 16-sided cells, 
respectively [21]. Kusner [22] has shown that a single cell with minimal interfaces 
in a froth which is locally Euclidean or hyperbolic cannot have less than 13.39... faces 
on average. It is also known that the minimal number of faces per cell of a periodic, 
monotiled froth is 14. Weaire and Phelan have recently given an example of froth with 
(/) = 13.5 (the so-called A15 phase) which minimizes the total interfacial area [18] (see 
also [23]). 

Natural froths minimize their free energy { Configurational Energymhms Temper- 
ature X Entropy). With the bounds given above, this condition is realized when the 
value of (/) is between 13.29... and 13.5 (or 14 for periodic monotiled froths). The lower 
bound corresponds to configurations with maximal orientational entropy, whereas the 
upper bound corresponds to configurations with minimal interfacial energy. 

There exists a class of natural structures, the Frank and Kasper phases (or the 
larger class of t.c.p. structures [16], [17]), for which (/) falls within these bounds. 
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(5.1) 



6- (n) 




(5.2) 



(/)* = 8 + 2v^ 
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These structures are periodic and made of 12-, 14-, 15- and 16-sided cells whose faces 
are either pentagons or hexagons. It can be verified that some of them fulfil the 
condition of Euclidean space-filling given by Eq.(4.1). We can therefore assume that 
the t.c.p. structures are Euclidean shell-structured- infiatable froths. Then, their shell- 
network is a periodic tiling made of pentagons and hexagons only. Let the 2D unit 
cell of the shell-network consist of /^^^ pentagons and /^^^ hexagons, belonging to N* 
polyhedra within the layer between two subsequent shells. The number of polyhedra 
in the 3D unit cell is a multiple of N*. The average number of edges per face in the 
shell-network is 

^^^^= /(6) + /(5) • (5.3) 
Substituting into Eq.(4.1), one obtains 

20/('^) + 14/(^) 

2/(6) • 

The number of polyhedra in the 3D unit cell can be calculated with the help of the 
numbers of faces the "outgoing" (/_|_) and "incoming" (/_) froths in the unit cell of 
the shell-network. These numbers coincide with the numbers of polyhedra in the layers 
above (/_|_) and below (/_) the shell which have one or more faces belonging to the 2D 
unit cell. In the limit of large shell-networks, one has the relation — 2/+(_), with 

vj^ (resp. V-) counting the number of 3-connected vertices in the 2D unit cell which 
belong to the "outgoing" (resp. "incoming") froths. Eq.(3.10) can then be written in 
terms of the quantities associated with the 2D unit cell only 

/, + /_=.jM^) (6.5) 



6- (n) 



N 



{vx counts the number of 4-connected vertices in the 2D unit cell). On the other hand, 
since the 3D system is periodic, one has /+-|-/- = 2A'"* on average. Therefore Eq.(5.5) 
can be written as 



V 



iV* = f(^l + 2^| . (5.6) 

If one puts into Eq.(5.4) the simplest combinations of integers /^^^ and J^^) which are 
such that (/) falls within the two bounds 13.29... and 13.5, one retrieves the average 
number of faces of the 3D unit cell of all experimentally known t.c.p., which are listed 
in Table 1. (The Table gives all the possible combinations {f^^\f^^^) up to Z^^) =4 
and, for /^^^ > 4, only those corresponding to known natural structures.) Also given 
are the corresponding values of N* , obtained from Eq.(5.6). These values of N* are 
exactly equal to the sum of the lowest non congruent numbers of 16- (p), 15- (q), 14- 
(r) and 12-sided polyhedra {x) in the structural formula of the corresponding t.c.p. 
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[17]. The Table presents also several simple combinations {f^^\ Z^^-*) which correspond 
to structures not (yet) observed (they are indicated by blanks in the last column). 
Notably, combinations (2,23), (2,25)... may be good candidates for t.c.p. structures 
yet to be observed. On the other hand, combinations (2,19), (3,28), (3,29) and (4,39) 
may be too distorted to qualify as t.c.p. structures. They may be realized with atoms of 
very different sizes. Note finally, that when (/) is represented as a function of the ratio 
y(5)^y(6)^ the structures in the Table 1 tend to gather into distinct groups. This may 
indicate either the existence of unfavourable configurations or structural mode-locking 
into the simplest t.c.p. structures (A15, Z, cr, ... C15). 

All these facts strongly suggests that the t.c.p. are shell-structured-inflatable 
froths. 

VI. CONCLUSION 

In this paper we have introduced a new way to study froths, which emphasizes 
their shell structure. We have studied an important subclass of shell-structured froths, 
i.e. those which can be generated in a recursive way according to an inflationary pro- 
cedure. For 2D froths (and networks with any coordination number) and 3D froths 
we have found that this recursive procedure is described by the logistic map. This 
map allows for a natural differentiation between froths tiling elliptic, hyperbolic or 
Euclidean manifolds, without any a-priori imposed curvature condition. In particular, 
the logistic map is able to characterize 3D curved manifolds, thereby providing a way 
to deflne the curvature from topological considerations when the Gauss-Bonnet theo- 
rem is not applicable. The logistic map in 3D enables us recover known space-fllling 
conflgurations, and also to suggest new ones. It is clear that the approach using the 
logistic map is very powerful, since classiflcation of 3D space-fllling conflgurations is 
reduced to the study of the 2D tiling of the (elliptic) shell-surface. As an example of 
the power and generality of this approach, we have been able to retrieve the topological 
properties of all experimentally-known t.c.p. structures by studying the tiling of the 
shell-surface by pentagons and hexagons. 
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APPENDICES 



A INFLATION OF TWO DIMENSIONAL Z-VALENT NETWORKS 
WITH Z > 4 

The generalization of Eq.(2.3) to the description of 2D shell-structured-inflatable 
froths with coordination number 2 > 4 is as follows. Every shell has {z — 1) different 
types of vertices. Extending the notation of Section 2, the various types of vertices are 
labelled by V^^^ with a = 0, 1, z-2, vj*^ being the number of vertices belonging to 
shell (t) from which a (resp. z — 2 — a) edges are pointing towards shell (t + 1) (resp. 

shell (t — 1)). Every vertex vi*^ adds a cells between shells (t) and (t + 1). The total 
number of cells F^*) between the two shells is 

FW = X:aV;(^) = X:(.-a-2)V;(*+^) . (Al) 

a=0 0=0 

Let (n) denote the average number of edges per cell layer (t). If one sums over all cells 
in this layer, one obtains 

z-2 z-2 

(n)FW = 5^(a + l)Vi*)+5^(^-a-l)yi^+i) . (A2) 



a=0 0=0 



Since 



and 



z — a 



one has 

z-2 



{A.5) 



a=0 a=0 
a=0 0=0 

One obtains finally the recursion relation 

((n) - 2) - 2(^ - 1)) = + . (A.6) 
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The matrix form of this recursion relation is the same as for z = 3 (Eq.(2.3) ), with 
recursion parameter s = {n){z—2)—2{z—l). The initial conditions are F^-*^) = {z—2){n) 
and = 1. 

Euclidean tilings are associated with the fixed point s* — 2, i.e. to the equation 

(n) = -^ . (A7) 
z — z 

The only regular solutions {z and (n) integers) of this equation are (6, 3) (tiling by 
triangles), (4,4) (tiling by squares) and (3,6) (tiling by hexagons) (dual of (6,3)). 

B.l SHELL-STRUCTURED BUT NON-INFLATABLE 2D FROTHS 

Some 2D shell-structured froths cannot be constructed according to the recursion 
procedure of Eq.(2.3). These froths have local inclusions which are topological defects 
in the recursion procedure. An inclusion in a layer is a cell with neighbouring cells in 
this layer and only in one of the two neighbouring layers. Topological defects fall in 
two classes: vertex decorations (Figs. (9a) and (9b)) and edge decorations (Fig. (9c)). 
In all cases the inclusion is on the + side of shell (t). 

Defects can be eliminated by removing one or more edges and its surrounding 
vertices. A vcrtcx-dccoration defect is then replaced by an ordinary vertex (Fig. (10a)). 
An edge-decoration defect is then replaced by edges on the shell (Fig. (10b)). 

The removal of one edge reduces by one unit the number of faces in the layer. 
This operation corresponds to the transformations E ^ E — 3, F— > V — 2, F— > F — 1. 
Consequently, since (n) = the average number of edges per cell changes as 

(n)'=(n) + -^((n)-6) . (S.1.1) 

The recursion parameter s = (n) — 4 changes therefore as 

s' = s + j^{s-2) . (S.1.2) 

One sees that the fixed point s* = 2 remains unchanged by defect elimination. More- 
over, elliptic froths become more elliptic (i.e. (n)' < (n) < 6) whereas hyperbolic 
froths become more hyperbolic (i.e. (n)' > (n) > 6). Thus, the Euclidean, hyperbolic 
or elliptic character of the manifold tiled by the froth is not modified by defect elim- 
ination (it is indeed given by the Euler-Poincare characteristic which is a topological 
invariant). 

B.2 SHELL-STRUCTURED BUT NON-INFLATABLE 3D FROTHS 

By analogy with the 2D case, one can define a topological distance r between two 
cells A and B as the minimal number of faces that must be crossed by a path that 
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connects A and B. A 3D shell-structured- inflatable froth is defined by the following 
two conditions : 

1) For any set of cells equidistant from a germ cell, there exists a closed non self- 
intersecting surface that cuts these cells and no others. 

2) Any cell at distance t from the germ cell is the neighbour of at least one cell at 
distance t + 1. 

Shells are closed surfaces tiled by faces of cells; they bound layers of equidistant cells. 
It is possible to connect two adjacent shells (t) and {t + 1) through a set of faces, each 
with one edge on shell (t) and one on shell (t+l). Shell (t) separates the whole froth 
into an inner froth, constituted of cells at distances r < t, and an outer froth, with 
cells at distances r > t. 

There are local defects which violate rules (1) or (2). These non-inflatable con- 
flgurations in 3D froths are shown in Fig. (11). These are particular examples of the 
three general classes of 3D topological defects: vertex, edge and face decoration, as in 
2D these non-inflatable configuration can be eliminated. Defects elimination is made 
by removing one (or more) facc(s), together with surrounding edges and vertices. The 
removal of one face with n edges reduces by one unit the total number C of cells. 
This operation corresponds to the transformation C — > C — 1 and F — > F — 1 — n. 
Consequently, since (/) = the average number of faces per cell changes as 

(/)' = (/) + ^ ((/) - 2(n + 1)) . (S.2.1) 

In contrast to the 2D case, this transformation depends on the parameter n. This is 
not surprising since it is well-known that in 3D the value of (/) is not directly related 
to the curvature of the manifold tiled by the froth. 

C RANDOM 3D EUCLIDEAN FROTHS FROM 
2D RANDOM SHELL NETWORKS 

Eq.(4.1) implies that a 3D random froth can be constructed from the superposition 
of two 2D random froths. To study this general case it is useful to rewrite Eq.(4.1) in 
term of the number of intersections of edges of the incoming froths by edges of the 
outgoing froths and vice versa. For a given shell t this quantity is equal to 

where we used the identity 3V'^*^''_^ = 2E^^^_y Using equation (3.10) it is possible to 
express in terms of {n)^- One has 



2 / 6-{n)N 4(n)jv \ 
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When the number of network cells is much larger than unity, one has = (2/3)(6 — 
{i^)N)/{{n)N — 4). Substituting into (4.1), one obtains 

(/) = 10 + . (as) 

In principle, in random froths, can take any value between zero and infinity 
(but only between 2/3 and 1 for periodic monotiled froths). For example, p^ = oo 
corresponds to a froth made with layers of infinitely long bricks disposed, layer by 
layer, with orientation alternating by 90°. In this case the network is a square lattice. 
The opposite limit {p^ = 0) corresponds, for example, to a 3D froth made with layers 
of large and small cells, when the ratio between the cell- sizes tends to infinity. In this 
case the network is the result of the superposition of a froth with cells of large sizes 
and a froth of small sizes and the probability of intersection of two edges of these two 
froths tends to zero. 
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FIGURE CAPTIONS 



Fig.l 

Schematic picture of a 2D shell-structured-inflatable froth. 
Fig.2 

Any 3D shell (t) is tiled with a network generated by the intersection of the faces 
coming to and going away from its surface. This shell-network has 4-connected vertices 
y^*^ (whose all 4 edges are belonging to the shell-network) and 3-connected vertices 
^+{-) (with 3 edges belonging to the shell-network and the last one going away from 
it)." 

Fig.3 

An example of 3D shell-structured-inflatable froth, the Kelvin froth. A portion of 
the shell-network is brought out by hatcheries. 

Fig.4 

The two 3D space-flUing unit cells constructed from the shell-network (a). The 
cell (b) is topologically equivalent to Kelvin's a-tetrakaidecahedron and the cell (c) to 
its twisted variant. Both have 14 faces. 

Fig.5 

The two 3D space-filling unit cells constructed from the shell-network (a) gener- 
ated by the superposition of two "squeezed" hexagonal lattices. The cell (b) is topolog- 
ically equivalent to the /3-tetrakaidecahedron. The cell (c) is topologically equivalent 
to the Goldberg cell. 

Fig. 6 

The 3D space-filling unit cell (b) (which has 16 faces) resulting from the shell- 
network (a) generated by the superposition of two "squeezed" hexagonal lattices. 

Fig.7 

Example of a 3D periodic shell-structured-infiatable froth (with (/) = 12) whose 
unit cell has two different elementary cells, (a) Shell-network, (b) 3D unit cell. 
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Fig.8 

The average number (/) of faces per cell in a froth plotted as a function of the 
average number {n)^ of edges per 2D cell in the shell network (Eq.(4.1), curve labelled 
"space-filling") and of the average number (n) of edges per face in the froth (Eq.(5.1), 
curve labelled "3D"). The abscissa n represents both {n)N and (n). 

Fig.9 

Local topological defects in the 2D recursion procedure. Figures (a) and (b) are 
examples of vertex decorations whereas Figure (c) is an example of edge decoration. 
The index t denotes the topological distance. 

Fig. 10 

Schematic representations of the elimination of a 2D local topological defect, (a) 
Vertex decoration, (b) Edge decoration. 

Fig. 11 

Local topological defects in the 3D recursion procedure, (a) is a vertex decoration 
defect, (b) is an edge decoration defect and (c) is a face decoration defect. The index 
t denotes the topological distance. 

TABLE CAPTION 

Table 1 

Average number of faces (/) and (minimal) number of elements in the 3D unit 
cell A^* of all the t.c.p. structures known experimentally (labelled in the last column) 
and of hypothetical t.c.p. structures (indicated by a blank in the last column). The 
integers p, q, r and x indicate respectively the proportions of 3D cells with 16, 15, 14 
and 12 faces present in the 3D unit cell. 
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Number of 
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